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The origin of radio haloes and non-thermal emission in 
clusters of galaxies 



O 
O 

(N 

< 
oo 



> 
OS 

o 
in 

00 

o 

O 

a: 
6 

c3 



H. Liang 1 , V. A. Dogiel 2 , and M. Birkinshaw 1 

1 Physics Dept., University of Bristol, Tyndall Avenue, Bristol BS8 1TL, England 

2 P. N. Lebedev Physical Institute, 117924 Moscow, Russia 



ABSTRACT 

We study the origin of the non-thermal emission from the intracluster medium, in- 
cluding the excess hard X-ray emission and cluster- wide radio haloes, through fitting 
two representative models to the Coma cluster. If the synchrotron emitting relativis- 
tic electrons are accelerated in situ from the vast pool of thermal electrons, then a 
quasi-stationary solution of the kinetic equation with particle acceleration through tur- 
bulence at high energies (> 200 keV) naturally produces a population of supra-thermal 
electrons responsible for the excess hard X-ray emission through bremsstrahlung. In- 
verse Compton scattering is negligible at hard X-ray energies in this case. The ra- 
dio halo flux density constrains the magnetic field strength to a value close to that 
of equipartition ~ 1/iG. Alternatively, if the relativistic electrons are injected from 
numerous localised 'external' sources then the hard X-rays are best explained by in- 
verse Compton scattering from GeV electrons, and little of the hard X-radiation has 
a bremsstrahlung origin. In this case, the magnetic field strength is constrained to 
~ 0.1 — 0.2 fiG. Both models assume that the non-thermal emissions are generated by 
a single electron spectrum, so that only two free parameters, well constrained by the 
observed hard X-ray and radio halo spectra, are needed in either case. Measurements 
of the cluster magnetic field will distinguish between the models. 

Key words: acceleration of particles - radiation mechanism: non-thermal - galaxies: 
clusters: general - radio continuum: general - X-rays: general 



1 INTRODUCTION 

Since the late 1960s, clusters of galaxies have been known to 
be strong X-ray sources where the thermal emission from the 
intracluster medium (ICM) produces a diffuse X-ray halo. 
The ICM is diffuse (central electron density of ~ 10 -3 cm -3 ) 
with typical temperature 2 x 10 7 to 2 x 10 s K. In the 1970s, 
similarly diffuse radio emission was found in the Coma clus- 
ter. This radio emission was found to have a steep power 
law index (a ~ 1 for S v oc v~ a ) which indicates a non- 
thermal origin. Despite efforts to detect similar radio haloes 
in other clusters, until recently few were found with diffuse 
radio emission (see review by Feretti & Giovannini 1996). 
This led to the belief that radio haloes are rare. However, 
the number of clusters found to have diffuse radio emission 
has more than doubled over the last couple of years, mainly 
because of the availability of moder ate brightness sensitivit y 
surveys at 1.4 GHz, e.g. the NVSS (|Giovannini et al. 1990|) , 



as well as high brightness sensitivity observations of hot clus- 
ters to detect the Sunyaev-ZePdovich effect (e.g. Herbig & 
Birkinshaw 1994; Liang et al. 2000). 

Diffuse cluster radio emission is now divided into two 
distinct classes: haloes and relics. Haloes permeate the en- 



tire cluster and resemble the thermal X-ray emission, but 
relics are on the cluster periphery. Both haloes and relics 
exhibit steep spectral indices (a ~ 1) and are extended 
over ~ 1 Mpc with the haloes having lower surface bright- 
ness. Emission from radio relics is strongly polarised (up to 
~ 20%) but no polarised emission has been detected in radio 
haloes. Radio haloes and relics are believed to emit by the 
synchrotron process. The non-detection of polarised emis- 
sion from radio haloes is attributed to low surface brightness 
and to Faraday depolarisation effects due to the mixing of 
thermal and relativistic plasma, especially towards the clus- 
ter centres. 

However, synchrotron emission from haloes requires a 
population of relativistic (GeV energy range) electrons and 
a cluster-wide magnetic field. The origins of both particles 
and field have been a subject of much debate since the late 
1970s (e.g. Jaffe 1977; Dennison 1980; Roland 1981). Re- 
cent evidence on this question has come through the ob- 
servation that radio haloes are preferentially found in high 
X-ray luminosity clusters (Giovannini et al. 1999), and that 
there may be a steep correlation between the radio power 
of haloes and the ICM temperature (e.g. Liang et al. 2000). 
No upper limits so far contradict this correlation in the tern- 



2 H. Liang et al. 



perature range 5-15 keV (Liang 2000). The relativistic par- 
ticles cannot simply leak out of radio sources in the clus- 
ter, because the diffusion path length (~ lOkpc assuming 
Alfven speed propagation) is short compared to the size of 
clusters (> IMpc). Giovannini & Feretti (2000) found that 
while diffuse radio sources were found in bright X-ray clus- 
ters, none of the 11 clusters selected for the presence of a 
tailed radio source within the central 300 kpc showed any 
diffuse radio emission. The presence of diffuse radio emis- 
sion is therefore more closely connected with the proper- 
ties of the thermal plasma in the ICM than the presence of 
tailed radio sources which were thought to provide relativis- 
tic particles. This gives observational evidence that particles 
escaping from radio sources in the cluster are unlikely to be 
accelerated to the relativistic energies required for the radio 
halo emission. Instead, we interpret the correlation between 
halo radio power and ICM temperature, and the similarity 
in appearance of the X-ray and radio haloes, as showing 
that the relativistic electrons may be accelerated from the 
thermal pool. The direct cause of this acceleration may be 
the merging of subclusters and turbulence in the ICM (e.g. 
Tribble 1993, Brunetti et al. 2001). 

The idea that the same relativistic electrons generate 
radio emission through their synchrotron losses and X-ray 
or gamma-ray emission by inverse Compton (IC) scattering 
was first used by Felten and Morrison (1966) for the inter- 
pretation of the galactic diffuse emission. Later, Cooke et 
al. (1978) used it for the interpretation of the X-ray emis- 
sion from the lobes of the radio galaxy Centaurus A, as- 
suming that the electrons are scattered by the cosmic mi- 
crowave background (CMB) photons. A similar model was 
used by Rephaeli (1979) for the interpretation of the hard 
X-ray emission from the Coma cluster: the relativistic elec- 
trons responsible for the radio halo emission can scatter the 
CMB photons up to X-ray energies through IC scattering. 
Efforts have been made to detect this excess of non-thermal 
emission above the thermal emission from the hot ICM. It 
is easiest to detect at hard X-ray energies (above 20keV) 
where the thermal spectrum is expected to cut off expo- 
nentially. Recent observations from Beppo-SAX have shown 
that such a hard X-ray excess does exist in some clusters, 
e.g. Coma and A2256 (Fusco-Femiano et al. 1999, 2000). 

However, alternative theories of the origin of a hard ex- 
cess have been proposed by various authors. Enfilin et al. 
(1999) first suggested that the hard X-ray excess observed 
in the Coma cluster may be the result of bremsstrahlung 
of supra-thermal electrons in the ICM accelerated by tur- 
bulence. Both Enfilin et al. (1999) and Sarazin & Kemp- 
ner (2000) have assumed the acceleration process works on 
particles with energies of a few times kT x (i.e. a few tens 
of keV), and that the electron energy distribution has a 
sharp transition from a Maxwellian to a power law at a 
few times kT x . Dogiel (2000) pointed out that the acceler- 
ation of particles is only efficient above ~ 900 keV (i.e. the 
power law distribution starts at ~ 900 keV), but the distor- 
tion of the Maxwellian is already significant at much lower 
energies (as low as ~ 30 keV) due to the flux of thermal 
particles running away into the region of acceleration (for 
details see Gurevich 1960). It is the 'quasi-thermal' particles 
just above ~ 30 keV that produce the observed hard X-ray 
excess through bremsstrahlung emission. However, our anal- 
ysis of the in situ acceleration model given below (Sec. 6) 



shows that the spectrum of MHD turbulence needed to pro- 
duce the radio emitting electrons is very steep which may 
make this model problematic. 

On the other hand, a model where the X-ray and radio 
emitting particles are generated by electrons from an 'ex- 
ternal' injected population cannot be ruled out. A possible 
model could be that particles in the intracluster medium are 
accelerated only in regions with high level of turbulence, for 
example near collisionless bow shocks generated by motions 
of hypothetical dark matter haloes (e.g., Bykov et al. 2000b). 
Radio observations have not shown significant spatial fluc- 
tuations of intensity in radio haloes (e.g. Giovannini et al. 
1993, Liang et al. 2000) which means that we must have ei- 
ther in situ acceleration of these electrons everywhere in the 
intracluster medium or injection by numerous sources whose 
separation is smaller than the path length of the electrons 
(~ 10 kpc) or the length corresponding to the resolution of 
the radio maps. In the case of the Coma cluster which is 
relatively bright and close by, the highest resolution VLA 
map of the radio halo corresponds to a length scale similar 
to the diffusion path length of the elecrons. Therefore, the 
injection model requires a high density of localised 'external' 
sources. 

Below we investigate both the local and global acceler- 
ation models and determine formally the model parameters 
required to fit the whole observed spectrum from radio to 
X-rays for the Coma cluster, since this is the best-studied 
cluster in all wave-bands. 

Our analysis is performed in the framework of a model 
where spatial variations of the model parameters are ne- 
glected. There are several reasons for this: a) observations 
show more or less uniform distribution of non-thermal emis- 
sion in the halo; and b) we do not have enough informa- 
tion about spatial variations of the hard X-ray flux from 
the Coma halo, therefore we derive averaged parameters of 
particle acceleration based on the total flux of hard X-rays. 

The cosmological parameters adopted are Ho = 
50kms _1 Mpc -1 , go = 0.5, so that for the Coma cluster 
1 arcmin is 40 kpc. 



2 THE COMA CLUSTER 

Extensive data from radio to X-ray bands exist for the Coma 
cluster (z ~ 0.0232). The best estimates of the ICM tem- 
perature are T x = 8.21 ± 0.09 keV from the Ginga satellite 
(Hughes et al. 1993) and T x = 8.25 ± 0.10 keV from the 
XMM-Newton satellite (Arnaud et al. 2001). The new high 
resolution XMM-Newton data have shown that at least the 
central 1 Mpc of the cluster is isothermal. The gas distri- 
bution is well parameterised by a /3-model (Cavaliere et al. 
1979) with electron density distribution given by 



n e (r) = n [l + (r/r a ) 2 ]- 3 ^ 2 

where /3 = 0.75 ± 0.03, the core radius ro 
408h~g kpc) 



(1) 

= 10.'5±0!6 (i.e. 
and the central 

3 J, 1 / 2 , 



electron density n = (2.89 ± 0.04) x 10~ 3 h 5 ' cm" d , from 
the ROSAT all-sky survey data (Briel et al. 1992). 

An excess of hard X-rays above the thermal component 
has been been detected in the 20-80 keV range by Beppo- 
SAX satellite ( Fusco-Femiano et al. 1999] ) . From Eq. []J we 



estimate the average electron density and squared electron 
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density within the Beppo-SAX field of view of 1?3 (i.e. a vol- 



ume of V ~ 1.7 x 10' b cm' to be n e 



1.23 x 10" 4 cm- d and 



nl ~ 5.22 x 10~ cm - . We adopt these values in the dis- 
cussions that follow, including the calculation of Fig. [I] & ^| 
Coma-C was the first radio halo detected, and its emis- 
sion has been measured between 30 MHz and 4850 MHz 
(Deiss et al. 1997 and references therein). The spectrum is 
consistent with a power law of index a ~ 1.25 between 30 
and 1400 MHz. The measurement at 4850 MHz is an upper 
limit and the apparent steepening at 2700 MHz has been 
contested (Deiss et al. 1997). Radio spectral index measure- 
ments for such diffuse sources are difficult and a mismatch 
of the beam size or difference in the brightness sensitivity 
of the observations can often lead to an apparent steepen- 
ing of the spectral index. Accordingly, we will disregard the 
measurements above 1400 MHz. The largest extent of the 
radio halo was found to be 80' x 45' on an Effelsberg map at 
1.4 GHz (Deiss et al. 1997); by contrast, the 1.4 GHz VLA 
image which is not as sensitive to the large scale emission 
showed a maximum extent of only 22' (Kim et al. 1990). 



3 ACCELERATION PROCESSES 

In most circumstances, the main source of charged parti- 
cle acceleration in a cosmic plasma is collisions of electrons 
and ions with magnetic field fluctuations. This leads to a 
slow, stochastic, energy gain. Here we consider magneto- 
hydrodynamic turbulence in the strong and weak limits. We 
can determine which regime is most appropriate in clusters 
from the spectral indices of radio haloes, since the spectral 
signatures of the two cases differ. 



3.1 Weak Turbulence 

If the magnetic turbulence is 'weak', i.e. the amplitude 
of magnetic field fluctuations, SB, is much less than the 
strength of the large scale magnetic field Bo, SB <C -Bo, 
then the interaction between the particles and the fluctua- 
tions has a resonance character, which has been discussed 
in the context of the interstellar medium (Berezinskii et al. 
1990). The resonance condition is 



uj(k) - =F = 0. 



(2) 



where the different signs correspond to the different wave 
polarisations, k and u are the wave-number and the wave 
frequency respectively, v» is the component of the parti- 
cle velocity along the magnetic field, and uih is the gyro- 
frequency of a particle with the total energy E to t and the 
charge Ze 



UJH = 



ZeB c 
Etot 



(3) 



The resonance condition (Eq. 0) can be written in the 



form 



ZeB 



pc(n — u)(k)/kv) 



(4) 



where v is the particle velocity, and /i is the cosine of the par- 
ticle pitch angle (i.e. the angle between the particle momen- 
tum, p, and the large scale magnetic field, Bo). For Alfven 
waves with dispersion relation 



uj(k) = ±\k\\\vA 



(•») 



where va is the Alfven velocity, and for fast particles (v 2> 
va) the resonance condition reduces to 



fen 



ZeB 
pen 



1 



r L \n\ 



(6) 



Here tl is the particle gyro-radius. We see that in general 
electrons and protons of the same energy (especially for non- 
relativistic energies) are scattered by waves with different 
wavelength and polarisation. Therefore, in some conditions 
electrons may be preferentially accelerated. This may be the 
case in the Coma halo; and in this respect models of hard X- 
ray emission from Coma which take into account of only the 
electron acceleration are not implausible (e.g. Schlickeiser et 
al. 1987). 

The process of acceleration may be described as mo- 
mentum diffusion with a diffusion coefficient determined by 
the spectrum of magneto-hydrodynamic turbulence. If the 
spectrum of waves is a power-law 



W(k) oc k' 



(T) 



the momentum diffusion coefficient depends on the particle 
momentum as (Berezinskii et al. 1990) 



a(j>) = a p q 



(8) 



It can be shown that in the simple thick target model 
the momentum distribution of accelerated particles is given 
by (see Gurevich 1960) 



fip) ocp 



-(9+1) 



(9) 



where f(p) is the electron number density distribution per 
momentum volume interval, so that the total number of par- 
ticles above a momentum p is given by 



N(>p) = / 4n f(p)p 2 dp, 



(10) 



and the value of the index q can be deduced from the spec- 
trum of non-thermal emission produced by the accelerated 
electrons. 

Below (Sec. 6) we perform the necessary analysis in an 
attempt to estimate ao and q and hence to define the mech- 
anism of particle acceleration in the intracluster medium of 
the Coma cluster. 



3.2 Strong Turbulence 

In the case of strong magneto-hydrodynamic turbulence, 
where the field fluctuations are of the same order as the 
mean field, both plasma components, positive and negative, 
are accelerated. This is seen, for example, in the accelera- 
tion of particles in molecular clouds (Dogiel et al. 1987) or 
acceleration by supersonic turbulence, which effectively acts 
as an ensemble of shock waves (Bykov and Toptygin 1993). 
Stochastic acceleration by strong turbulence is described as 
momentum diffusion with diffusion coefficient a(p) of the 
form (e.g. Toptygin 1985) 



a(p) = a p 



(11) 



which creates a relatively flat spectrum of accelerated par- 
ticles 
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fip) « p 



(12) 



The interpretation of the hard X-ray emission in Coma in 
terms of this kind of acceleration was discussed by Dogiel 
(2000). 

A similar acceleration process takes place near shock 
fronts, where the magnetic field strength jumps by an 
amount of the order of the mean field. However, the spec- 
trum of accelerated particles near a shock front has 

Hp) « p" 4 , (13) 

steeper than that of an ensemble of shocks (Eq. This 
case has been analysed for the Coma cluster by Bykov et al. 
(2000b). 



4 PARTICLE KINETIC EQUATION 

To describe the distribution function in both non-relativistic 
and relativistic energy ranges it is more convenient to use 
a canonical momentum variable, p, than the particle kinetic 
energy E. The electron number density distribution per unit 
energy interval, F(E), is related to the momentum distribu- 
tion per unit momentum volume, f(p), by 



F{E) = ^f{p)p 2 ^. 
where 



dp _ 1 

dE ~ p^m e c 2 kT x 



p 2 + 



(14) 



(15) 



where m e is the electron rest mass, T x is the temperature 
of the thermal gas, k is the Boltzman constant, and p is the 
dimensionless particle momentum, p — p j \J m<MT x where p 
is the dimensional momentum. 

The generalised Fokker-Planck approximation to the ki- 
netic equation, describing the particle distribution function 
/ in the intracluster medium (see Berezinskii et al. 1990), 
written in dimensionless variables is 



Of 



dt 
p 2 dp 



+ u(r) 



A(p) 



df _ p du(r) df 
dr 3 

df(p) 



dp 



dr dp 

+ B{p)f 



■Q(p,r) = 



■V(«(r,p)V/). 



(16) 



where u(r) is the particle convection, k is the spatial diffu- 
sion coefficient, and Q describes the distribution of sources 
emitting fast particles into the intracluster medium. We de- 
fine non-dimensional variables t = t • vo, r = r/L, where 
t and f are the dimensional variables, L is the size of the 
cluster, and we define the characteristic frequency i/ as 



vo = 



\fm e kT x 



(17) 



n e is the average electron number density, and r e is the 
classical electron radius r e = e /m e c . For the Coma halo 
of n e = 1.23 x 10~ 4 cm~ 3 within a volume of 1.7 x 10 75 cm -3 
and T x = 8.2 keV (Sec. 2), v ~ 6.4 x 10~ 17 s _1 . 

The frequency vo is proportional to the frequency v\ of 
thermal particle Coulomb collisions at kinetic energy kT x 



2v 



rn e c 
kT x 



In A, 



(18) 



where In A is the Coulomb logarithm. The frequency v\ ~ 
5.3 x 10~ 14 s" 1 for the Coma halo. 

In Eq. [n| B describes the particle energy gain and loss 
and A the momentum diffusion due to Coulomb collisions 
with background particles and interactions with the mag- 
netic fluctuations, 



B{p) = p 2 



A(p)= P 2 



(dp\ _ fdp\ _ (dp\ 

\dt J i on \dt J synIC V dt J 



'dp 



dp^ 



(-) 

V dt / ion 



Vt 5 



kT x 
m e c 2 



(19) 



(20) 



where 7 = 1 + E(p)/m e c 2 is the Lorentz factor, and the 
dimensionless parameter ag is related to the dimensional 
parameter fio as 



d {^Jm e kT x ) q - 2 



(21) 



As one can see from Eq. |2C)| the process of momentum 
diffusion is determined by Coulomb collisions with back- 
ground particles (the first rhs term of Eq. ^) and collisions 
with electromagnetic fluctuations (the second rhs term of 
Eq|H). 

The dp/dt terms are rates of momentum loss due 
to ionisation, synchrotron plus inverse Compton and 
bremsstrahlung emission respectively (e.g. Hayakawa 1969). 
The ionisation loss term is 



-_i / 2 | mec2 

\dtJion p\ kT x yf. 



7 



lot 



/ E(p)m e c 2 {-y 2 



\ h 2 ujp^ 2 



r-i 

0.43 



(22) 



where the plasma frequency lo p — \J 4-Ke 2 n e /m e . The com- 
bined synchrotron and IC loss term is 



V dt ) synl 



m e c* 



32-Kri 



1 

p\/ P + 'kT x 9Vm e kT x vo 
\m e c 2 J 



(23) 



where U cm b is the energy density of the CMB, and U mag is 
the magnetic energy density. The bremsstrahlung loss term 
is given by 



(I) =~v£+ 

V dt J brem p V 



2 - 2 2 2 i 

m e c z n e m e c e r e 1 
kT x hc\/m e kT x Vo 



(24) 



lof 



2(E(p) + m e c 2 ) 



E(p) 



In comparison to the loss terms, the average rate of acceler- 
ation due to momentum diffusion (Eq. 0) is 



\dtj c 



aop 



q-l 



(25) 



These loss terms and the acceleration term are plotted 
as a function of particle energy in Fig. |l] for the Coma clus- 
ter. At low energies ionisation loss dominates, and at high 
energies either the acceleration term or the combined syn- 
chrotron and IC term dominates depending on the value of 

q- 
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Figure 1. The average rate of momentum loss and acceleration as a function of energy in the Coma cluster (n e ~ 1.23 X 10 4 cm 3 and 
n\ ~ 5.22 X 10~ 8 cnr 6 within the Beppo-SAX field of view). The solid curve gives the rate of acceleration, where ao = 0.43 and q = 4.5 
corresponding to the best fit to observational data (see Sec. 6); the dotted curve gives the rate of ionisation loss; the dashed curve gives 
the rate of synchrotron and inverse-Compton loss; the long-dashed curve gives the rate of bremsstrahlung loss. For comparison with 
Dogiel (2000), we also show the acceleration rate corresponding to strong turbulence with ao = 5 and q = 2. 



Equation |l^ describes the particle spectrum in the en- 
ergy range from thermal energies, where the spectrum is 
Maxwellian, up to non-thermal energies, as well as the par- 
ticle spatial distribution. In general it is difficult to obtain a 
solution, though in some special cases it is possible to solve 
it analytically or numerically (see Bulanov and Dogiel 1979). 

Fortunately, in many cases the spatial distribution is 
unimportant. Whether or not a model without spatial terms 
can be applied to Coma is rather questionable. At the mo- 
ment we do not have enough information to prove or dis- 
prove this. 

If we ignore the spatial variation, then we can integrate 
Eq. [l(] over the volume of the object and the terms de- 
scribing the spatial distribution / converge to a single term, 
which can be represented as //t e ff where t e ff is the effective 
escape time of the particles due to their spatial propaga- 
tion. If the characteristic escape time is much greater than 
the energy loss timescale, we can neglect the term f/t e «- 
Then Eq. |l(| rewritten for the spatially-averaged distribu- 
tion function f(p), can be presented as 



dt 



1__9 

p 2 dp 



A( P ) 



am 

dp 



B(p)f 



= Q(p) ■ 



(26) 



Two possible fast particle production models are anal- 
ysed for the Coma halo in Sec. 6 and 7: 

• in situ acceleration everywhere in the halo, where the 
spectrum of non-thermal particles is established by accel- 



eration processes represented in Eq. ^ by the momentum 
diffusion term A(p) and balanced by the loss term, B(p); 
and 

• injection into the intracluster medium from many rel- 
atively small regions where acceleration takes place (per- 
haps strong shocks, or regions with high levels of magnetic 
turbulence). These regions are supposed to be more or less 
uniformly distributed in space with no acceleration between 
them. The spectrum in the intracluster medium is then fully 
determined by energy loss processes and the injection spec- 
trum, i.e. B(p) and Q(p)- 



5 PROCESSES OF NON-THERMAL EMISSION 
GENERATED BY FAST ELECTRONS 

To determine the parameters for both models we need to 
calculate the X-ray and radio spectra expected from the 
non-thermal electrons and compare our results with obser- 
vational data. In this section, we collect the equations de- 
scribing the radiation processes. 



5.1 Non-thermal Bremsstrahlung Emission 

The bremsstrahlung flux density Ft rem in units of 
ph cm -2 s _1 keV -1 is given by 
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v 



JE _ 



m e dE-y 



(27) 



where V is the emission volume, cLl is the luminosity dis- 
tance, is the photon energy (not to be confused with the 
particle energy E), and the emission cross section is given 
by 



da x 16 e 2 r 2 e m e c 2 , ( VE + y/E-E^ 
■ log ' 



(28) 



(Hayakawa 1969). Non-thermal bremsstrahlung radiation in 
the hard X-ray range comes from electrons in the energy 
range ~ 30 - 100 keV. 



5.2 Inverse Compton Emission 

Here we calculate the IC flux density Fic in units of 
ph cm -2 s _1 keV -1 due to the scattering of CMB pho- 
tons by relativistic particles. Since we are only interested 
in the IC scattering above 100 eV (i.e. EUV and above), we 
are interested only in electrons with E > 10 s keV. We will 
show (Fig. ^) that the electron energy distribution follows a 
power law with slope 1 — q for E > 10 4 keV. We can use the 
analytic solution for a power law distribution in the Thom- 
son limit (a/ E^kT r <C m e c 2 ) given in Blumenthal & Gould 
(1970; Eq. 2.65) for this power law to find 



Fic 



V 



4nd1 7r/i 3 c 2 



K e {m e (?) 2 q 



(kT r ) q -^C(q- 



l)E~ q/2 ,(29) 



where numerical values of (^(q — 1) are given in Table 1 of 
Blumenthal & Gould (1970), T r is the temperature of the 
CMB (not to be confused with the electron temperature 
T x ), K e = F(E min )E mi „/(q — 2) is the normalisation for 
the power law distribution and E m in is typically 10 4 keV. 
Electrons of energies ~ 5 GeV inverse Compton scatter pho- 
tons of the CMB to hard X-ray energies ~ 60keV. 



5.3 Synchrotron Radio Emission 

We calculate the synchrotron flux density in Jy from the 
electron energy distribution using Ginzburg & Syrovatskii 
(1964): 

y rE max 

FraMo = — ?rl0 23 / V3r e eBF(E)f x dE, (30) 



where 



/• \ J sin 2 / A - -, ,0/)<VW- 



(31) 



v c = 4.21B sin6(E /m e c 2 ) 2 MHz, and B is in units of Gauss. 
Since in our case the electron energy distribution is a power 
law for E > 10 4 keV and the synchrotron emission in the 
radio frequency range is given by the GeV electrons, we can 
safely use the analytic equation for a power law electron en- 
ergy distribution given in Eq. 6.10 of Ginzburg & Syrovatskii 
(1964). For radio emission above ~ 10 MHz, we need elec- 
trons above energies of ~ 1 GeV for B ~ 1 fiG. 



6 MODEL OF IN SITU ACCELERATION 

If the injection term, Q(p) = 0, then apart from the ioni- 
sation and radiation losses, the spectrum is completely de- 
termined by two processes: Coulomb collisions, which form 
the thermal and quasi-thermal part of the spectrum, and 
stochastic acceleration, which forms the spectrum at high 
energies (see Fig. [jj). 

A quasi-stationary solution to Eq. ^6|, in the sense that 
the presence of a 'weak' acceleration mechanism leads to the 
continuous acceleration of particles above a certain injection 
energy, is given by Gurevich (1960): 



o A(v) 



dv g 



where 



9 = 



(32) 



(33) 



The energy distribution of the electron number density 
thus calculated is shown in Fig. ^ for 3 sets of (ao,q). It can 
be seen that F(E) starts to deviate from a pure Maxwellian 
at a few tens of keV and reaches a power law at ~ 1—10 MeV. 

Radio spectra implied by these (ao,q) parameter sets are 
shown in Fig. ^, and the hard X-ray spectra are shown in 
Fig. From Eq. it follows that the momentum spec- 
tral index of the accelerated particles is q — 1, hence q 
is related to the the spectral index a of the synchrotron 
emission as 9 = 2a + 2. The observed radio spectrum 
(Fig. |^) constrains the electron momentum spectral index 
to q ~ 4.5lg'2 5 . Figure [5] shows that the best fit acceler- 
ation parameters from the radio and hard X-ray data are 
ao — 0.43^018, q — 4.5 =F 0.2, where the errors correspond to 
95% confidence. The corresponding magnetic field derived 
from adjusting the model flux density to fit the data at 
326 MHz, is B ~ 1.51±^|^ fj,G. The magnetic field is well 
constrained in this model (with only 2 free parameters) by 
the X-ray and radio data, and is consistent with the limit de- 
duced from Faraday rotation measurements, B < 5.9/igg 2 /xG 
on field tangling scales of < lkpc (Feretti et al. 1995). In 
comparison, the equipartition magnetic field is B ~ 1 /iG 
assuming that there are equal number of relativistic protons 
and electrons and that the emission filling factor is 1. Any 
independent estimates of the cluster global magnetic field 
would provide a test for the validity of the model. 

In the framework of this model the origin of the hard 
non-thermal X-rays is bremsstrahlung since the flux of IC X- 
ray photons is negligible. Figure Q shows the expected ratio 
of the bremsstrahlung to IC flux density for the Coma cluster 
at 70 keV as a function of the momentum diffusion coefficient 
spectral index, q. If the GeV electrons responsible for the ra- 
dio halo emission are the result of in-situ acceleration, then 
the observed radio halo spectrum (a ~ 1.25 or q ~ 4.5) 
implies that the contribution from the bremsstrahlung of 
30-100 keV electrons is always greater than that of the IC 
scattering of CMB photons by the GeV electrons. This re- 
sult is independent of the magnetic field, since neither the 
bremsstrahlung nor the IC emissivity depends on the mag- 
netic field. 

CMB photons can also be up-scattered to the extreme- 
ultraviolet (EUV) energy range by electrons of energy 
2.5 x 10 J keV through inverse Compton scattering. The 
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Figure 2. Energy distribution of electron density for a pure thermal distribution (solid curve), and various acceleration parameters (ao) 
and spectral indices q. The dotted line is for ao = 0.3&g = 5; the dashed line is for ao = 0.43 & q = 4.5; and the long-dashed line is 



for ao = 0.6 & q = 4. The cluster is assumed to be isothermal with T x = 8.2 keV and average electron density n e 



1.23 X 10- 



within the Beppo-SAX field of view. The electrons of order ~ 10 7 keV are responsible for the radio emission and IC emission in the hard 
X-ray regime, and the electrons at ~ 50 keV give rise to non-thermal Bremsstrahlung emission in the hard X-rays. 



detection of excess EUV emission above that of thermal 
bremsstrahlung has been reported in the Coma cluster using 
data gathered by the Extreme Ultraviolet Explorer (EUVE; 
Lieu et al. 1999 and references therein). However, this re- 
sult is controversial (see e.g. Bowyer et al. 1999), hence we 
have chosen to predict the excess EUV flux density in the 
framework of our model instead of fitting it to the data. We 
estimate the expected excess EUV flux density at 100 eV to 
be 1.9 x 10~ 3 ph cm 2 s 1 keV 1 within a 15 radius which 
is nearly 2 orders of magnitude below the EUVE sensitivity. 
However, given the controversy over the observational data, 
it can not be used to rule out the model without further 
independent measurements. 

It is interesting to note that, for q < 4, IC starts to dom- 
inate over bremsstrahlung in the hard X-ray regime, but for 
q < 2 bremsstrahlung dominates again because synchrotron 
and IC losses exceed the gains from acceleration and causes 
the electron spectrum to cut off at ~ 1 GeV. To produce 
X-ray emission above 20 keV through IC scattering of the 
CMB, we need electrons with energies above ~ 3.5 GeV. For 
example, if the radiating particles are accelerated by strong 
turbulence, that is q — 2 (hence ao = 5 to fit the hard X- 
ray data assuming a bremsstrahlung origin) , then ionisation 
losses dominate over acceleration below 1 MeV (see Fig. []]). 
The gain from acceleration overwhelms particle energy losses 
in the energy range of 1 MeV - 1 GeV. Hence, the maximum 
energy of accelerated particles is ~ 1 GeV. At higher en- 



ergies, the combined synchrotron and IC losses dominates 
over acceleration (see Fig. ^), hence strong turbulence is un- 
able to produce radio emitting GeV electrons. Acceleration 
by strong turbulence can produce only relatively low energy 
electrons emitting in the energy range of hard X-rays to 
soft gamma-rays through bremsstrahlung. Another mecha- 
nism for the production of radio-emitting electrons is then 
necessary to explain the radio halo in Coma. This conclu- 
sion is somewhat in contradiction with the results of Dogiel 
(2000) who showed that in this case strong turbulence can 
accelerate particles up to energies of ~ 10 GeV, albeit with 
a much flatter spectrum than needed to produce the ob- 
served radio spectrum. The reason for this disagreement is 
that Dogiel (2000) used the density for central regions of the 
cluster, where the frequency of thermal particle collisions is 
high. This led to a higher acceleration efficiency than in our 
present treatment which uses average cluster parameters. 

It is shown in Fig. ^ that for q < 2, the synchrotron 
and IC losses win over the acceleration term at energies 
even lower than 1 GeV. This explains why Blasi (2000), in 
his study of stochastic acceleration in the ICM assuming a 
Kolmogorov spectrum (q — 5/3) concluded that in situ ac- 
celeration could not explain the radio emission because of 
the synchrotron and IC loss cutoff. However, we have shown 
that the value q = 4.5, corresponding to weak turbulence 
(Sec. 3.2), is consistent with the radio data (Fig. ^). 

In the case of q > 3, however, the acceleration term 
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Figure 3. Synchrotron flux density derived from the electron energy distributions in Fig. ^ for the 3 sets of acceleration parameters 
(ao,q) and their corresponding magnetic field B. The magnetic field was adjusted such that the model flux densities fit the data at 
326 MHz. The observed flux density for the Coma radio halo are shown along with the 3 models. The solid line corresponds to q = 4.5 
(B = 1.51/iG); the dotted line corresponds to q = 5 (B = 9.7/xG); and the dashed line is for q = 4 (B = 0.16/xG). The average density 
within the radio emitting region (i.e. an equivalent radius of ~ 34') is ~ 1.4 X 10 — 4 cm - 3 . 



always wins over the loss terms at high energies (Fig. [lj), 
which means a rapid decrease of the characteristic accelera- 
tion time, Tacc, with particle momentum 



a(p) 



a u 



(34) 



and a larger q results in a smaller value of the injection en- 
ergy at which processes of acceleration dominate. The best 
fit parameters, ao ~ 0.43 and q ~ 4.5, imply a cluster aver- 
age T a cc ~ 7 x 10 6 yrs at the injection energy (~ 200 keV), 
which means that for E > 200 keV, r acc is much less than 
the cluster lifetime. Hence, the model is self-consistent in the 
sense that a 'quasi-stationary' state can easily be reached 
within a cluster lifetime. 

We note, however, the steep spectrum of turbulence 
(g ~ 4.5) required for the radio data may pose a problem 
for the model, since it is known that turbulence with q > 2 
cannot effectively scatter particles at pitch angles near 7r/2 
in the quasi- linear regime. Nonlinear processes have to be 
involved to cause effective acceleration(e.g., Berezinskii et 
al. 1990). Further consideration of this point is beyond the 
scope of our model. 

Cluster thermal X-ray and radio structures are similar 
in morphology (Deiss et al. 1997, Feretti 1999, Liang et al. 
2000, Govoni et al. 2001). This would be expected from the 
in situ model. We can calculate the expected relationship be- 
tween the radio halo and X-ray profile, if the synchrotron- 



emitting relativistic particles are drawn from the thermal 
pool through in-situ acceleration. If the magnetic field is 
produced by a dynamo process (e.g. turbulence), then mag- 
netic energy is in equipartition with the energy of turbulent 
motions (Kraichnan & Nagarajan 1967, De Young 1980). 
Assuming that the energy of turbulence is proportional to 
the thermal energy, then we have the magnetic energy den- 
sity in a cluster medium proportional to the thermal energy 
density P ma g oc P t h, so that 



B oc nl /2 T^ 2 . 



(35) 



As pointed out by Enfilin et al. (1998), this implies that 
the radio halo surface brightness profile should be slightly 
steeper than the X-ray profile. The X-ray surface brightness 
is given by the projection of n 2 ., which for a /3— model with 
n e given by Eq. |l| is 



Ix(r) oc 



l + (-) 2 

ro 



-3/3+1/2 



(36) 



for fixed T x . 

Since relativistic particles are accelerated from the ther- 
mal pool, we can assume that their density is proportional to 
the thermal particle density n e . Then the radio halo bright- 
ness profile is given by 



/l.4GHz oc / n e B 



dl. 



(37) 
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Figure 4. Models for hard X-ray (10-100 keV) emission in the Coma cluster. The solid curve shows a pure thermal bremsstrahlung 
spectrum calculated from the Raymond-Smith code (Raymond & Smith 1977) for T x = 8.2 keV; the dashed curve corresponds to the 
non-thermal bremsstrahlung emission derived from one of the electron energy distributions shown in Fig. ^ (dashed line; ag = 0.43 & 
q = 4.5); the dotted curve shows the X-ray emission from inverse-Compton scattering of the CMB by the rclativistic electrons of the 
same electron distribution. The data points correspond to the observed hard X-ray flux densities from Beppo-SAX. 



when using Eq. 
form of Eq. [l] 



and assuming an isothermal ICM of the 



AiGHz OC 



l + (-) 

ro 



-3/3(3+a)/4+l/2 



(38) 



If the radio spectral index is exactly a = 1, then the 
radio and X-ray profiles follow each other. In the case of 
Coma and 1E0657-56, the spectral index is a ~ 1.25 and the 
radio profile is slightly steeper than the X-ray profile. For 
the clusters studied by Govoni et al. (2001), the radio and 
X-ray surface brightness profiles are very similar in shape, 
implying a ~ 1. 



7 ENERGY OUTPUT OF NON-THERMAL 
AND QUASI-THERMAL PARTICLES 

In this section, for clarity we need to distinguish particles 
from the 3 different energy ranges: thermal, quasi-thermal 
and non-thermal. The Maxwellian thermal particle spec- 
trum is shaped by Coulomb collisions, where each thermal 
particle collides with other background particles, exchanging 
energies without changing the shape of the energy spectrum. 
The lifetime of the thermal particle distribution is infinite if 
we neglect thermal bremsstrahlung emission. On the other 
hand, a non-thermal particle spectrum is formed by accel- 
eration processes. If these particles are injected into a re- 
gion with no further acceleration, the particles would lose 



their energy to the background particles through Coulomb 
collisions, thus heating the background plasma and become 
thermal particles themselves. The lifetime of non-thermal 
particles is equal to that of the individual particles. Unlike 
the thermal and non-thermal particles, the quasi-thermal 
particles exist only in regions of active particle acceleration. 
Like the thermal particles, the quasi-thermal (or suprather- 
mal) particle distribution is formed by Coulomb collisions, 
but unlike the thermal particle distribution it is not in 
strict equilibrium (but in quasi-equilibrium) and the quasi- 
thermal particle distribution would dissipate like the non- 
thermal particles once the acceleration process is switched 
off. Coulomb collisions cause a net run-away of thermal par- 
ticles into regions of acceleration in the spectrum, resulting 
in a spectral distortion with an excess above the thermal 
equilibrium distribution in the quasi-thermal energy range. 
Hence the life-time of quasi-thermal particles is shorter than 
the thermal particles, but longer than the non-thermal par- 
ticles. 



The energy output of emitting particles is determined 
by their energy losses: ionisation losses (Coulomb collisions 
with background particles) for quasi-thermal particles, and 
inverse Compton scattering and synchrotron radiation losses 
for non-thermal electrons. The necessary energy output, L cr , 
can roughly be estimated from the observed value of X-ray 
luminosity L x as 
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Figure 5. The acceleration parameter ao and particle spectral index q for Coma based on the Beppo-SAX hard X-ray data and radio 
halo spectral index. The permitted ranges are ao ~ 0.25 — 0.7 and q ~ 4.3 — 4.7 with best fit values ao ~ 0.43 and q ~ 4.5. The dotted 
curve gives the best fit values of (ao,q) based on hard X-ray data, and the solid curves give the 95% confidence range. Similarly, the 
dashed line gives the best-fit q value based on the radio halo spectral index, with the solid lines giving the upper and lower 95% bounds. 



T 1±- 



(39) 



where t x is the characteristic time of X-ray emission by par- 
ent particles and r cr is the life-time of these particles. 

If Coma cluster emits X-rays by inverse Compton pro- 
cess, then from Eq. || it follows that 

4 ■ 10 43 erg s" 1 , (40) 



- tj ' r.r 



where L x ~ 4-10 erg s~ is the luminosity of hard 
(20 — 80 keV) X-rays from the Coma halo. The uncertainty 
in the estimate L^ r is determined by the unknown rate of 
synchrotron losses, since the magnetic field strength is only 
known within fairly wide bounds (a range that spans nearly 
2 orders of magnitude). 

In the case of a bremsstrahlung origin of the Coma 
luminosity produced by non-thermal particles, as in mod- 
els of Enfilin et al. (1999) and Sarazin & Kempner (2000), 
the necessary energy output L^ r is much larger than L x . 
The reason is that these non-thermal particles lose their en- 
ergy by collisions with background electrons and protons 
(ionisation losses) and the probability of their generating a 
bremsstrahlung photon is very low (< 2.5 x 10~ 4 for photons 
with E x < 100 keV; Valinia and Marshall 1998). Therefore 
the necessary energy output in this case is huge 

L° ~ 10 48 erg s" 1 , (41) 

and as pointed out by Petrosian (2001), this makes the 
model problematic. 

Similar problems arise in the interpretation of the hard 



X-ray flux from the galactic ridge. The spectrum of hard 
X-ray emission from the galactic ridge also shows an ex- 
cess above the thermal flux. However, in this case its origin 
cannot be interpreted as due to inverse Compton of rela- 
tivistic electrons (Skibo et al. 1996), thermal emission of 
hot plasma or the integrated contribution of discrete un- 
resolved source (Tanaka 2002). The only possibility to ex- 
plain the origin of hard X-ray flux from the galactic disk 
is just bremsstrahlung radiation of sub-relativistic galactic 
electrons. If these electrons are non-thermal the necessary 
energy output should be of the order of 10 43 erg s _1 , i.e. 
more than can be supplied by supernovae, the most proba- 
ble sources of energy in the Galaxy (see e.g. Valinia et al. 
2000, and Dogiel et al. 2002a). It was shown, however by Do- 
giel et al. (2002b) that this energy problem can be resolved 
if the emission is produced by quasi-thermal electrons. The 
life-time of quasi-thermal electrons is longer than that of 
non-thermal electrons. Therefore, as follows from Eq. |3s], the 
necessary energy output of quasi-thermal electrons is smaller 
than the energy output of non-thermal electrons. We showed 
above that in the case of in-situ acceleration, unlike the mod- 
els of Enfilin et al. (1999) and Sarazin & Kempner (2000) 
where non-thermal electrons with power-law spectrum were 
responsible for the excess hard X-ray in the Coma cluster, 
it is the quasi-thermal electrons with exponential distribu- 
tion that give rise to the excess hard X-ray emission. If our 
model is correct, the energy output in quasi-thermal elec- 
trons needed to explain the Coma X-ray flux is of the order 
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of Lg ~ 10 45 erg s 1 as follows from numerical calculations 
for the parameters ao = 0.43 and q — 4.5, i.e. much less than 
the estimate in Eq. |^ 

The bremsstrahlung X-ray flux can also be generated 
by sub-relativistic protons (see Dogiel 2001). Though we do 
not know which specific process is responsible for particle 
acceleration in Coma, we note that in the case of Fermi 
acceleration the characteristic acceleration frequency is pro- 
portional to ao ~ v/L (see Toptigyn (1985)) where v is the 
particle velocity and L is the particle mean free path (an av- 
erage distance between scattering centres). Since velocities 
of thermal electrons are yj ' M/m times larger than velocities 
of thermal protons, we expect that the rate of acceleration 
for electrons a§ > ajj, hence background electrons may be 
accelerated to sub-relativistic energies more effectively than 
protons. 



the tangling scale of the magnetic field in the Coma clus- 
ter, using an extended polarised radio source in the centre 
of the cluster, suggested that the cluster magnetic field in 
Coma has a small scale (< lkpc) component of strength 
<5.9/i5Q 2 /iG as well as a uniform large scale (~ 200 kpc) 
component of strength 0.1 — 0.3 pG which is consistent with 
the magnetic field of ~ 0.15 /iG deduced from the IC model. 

Figure ^ shows that for this model IC emission domi- 
nates bremsstrahlung at 70 keV for the required steep spec- 
tral index, rj ~ 4.5. The presence of strong shocks in galaxy 
clusters shown by spatial and temperature images from 
ROSAT, ASCA and now XMM-Newton and Chandra (see 
Donnelly et al. 1999, Markevitch et al. 1999, Arnaud et al. 
2001, Sun et al. 2001 and Boute 2001) lends observational 
support to this model. 



8 MODEL OF 'EXTERNAL' SOURCES OF 
ACCELERATION 

As follows from Fig. |l| the particle spectrum is controlled by 
ionisation losses at low energies and synchrotron+iC losses 
at high energies, which balance the source injection spec- 
trum and acceleration. If we assume no in situ acceleration, 
i.e. A(p) = in Eq. then the kinetic equation for the 
distribution function / can be written as: 



The stationary solution of this equation is given by 



f(p) = 



f(p) 



B( P ) 



P Q{p)dp 



K 



b{ p ) \(v-iy 



(42) 



(43) 



(44) 



where we assume a power law injection spectrum Q(p) = 
K"p~\ or Q(E) oc E~~ m (where 70 = r\ — 2 as follows from 
Eq.0). 

The best fit to the radio spectrum gives r\ ~ 4.5, or 
70 ~ 2.5 (Figure ^), just as is expected for the model of par- 
ticle acceleration by strong shocks (Biermann et al. 1993) . In 
this model, the origin of the hard X-ray emission is definitely 
IC with no significant bremsstrahlung (Fig. Q). The model 
reproduces the emission spectrum from hard X-rays to radio 
if the strength of the magnetic field is less than 0.1 — 0.2 fiG. 
Thus the finding from Faraday rotation measurements (e.g. 
Feretti et al. 1995; Clarke et al. 2001) of cluster fields of a 
few fj,G argues against the IC model for the hard X-ray (e.g. 
Enfilin et al. 1999 and Dogiel 2000). However, this is not 
necessarily a problem. First, since there are only a handful 
of radio sources behind the cluster with sufficient polarised 
flux to be detected, the regions of a cluster probed by these 
sources may not be representative of the average field of the 
entire cluster. Second, since as many as half of the polarised 
radio sources used to deduce the cluster magnetic field are 
cluster members rather than background sources, we expect 
that such measurements would be biased, as a higher B-field 
is likely near a radio source because of its perturbing effect 
on the ICM. Third, Feretti et al. (1995) in their study of 



9 CONCLUSIONS 

We have shown that it is possible to explain both the excess 
hard X-ray data and the radio halo spectrum with a model 
where the relativistic electrons are weakly accelerated from 
the thermal pool by turbulence in the ICM. By using a quasi- 
stationary solution of the kinetic equation, we arrive at an 
electron energy distribution that deviates from a Maxwellian 
at ~ 30keV and reaches a power law at ~ 1 — lOMeV. 
In the process of reaching an equilibrium, the collisions be- 
tween thermal particles produce the supra-thermal electrons 
at ~ 30 — 100 keV required to give the excess hard X-ray 
emission through bremsstrahlung emission. The accelerated 
GeV electrons are then sufficient to produce the observed ra- 
dio halo spectrum with a magnetic field close to the equipar- 
tition value. The model produces a self-consistent electron 
energy distribution and has only two free parameters, the ac- 
celeration parameter ao and the momentum spectral index 
q, which are well constrained by the observed excess hard X- 
ray spectrum and the radio halo spectral index (Fig. ^). In 
the hard X-ray energy range, the contribution from the in- 
verse Compton scattering of the GeV electrons is negligible 
compared with the non-thermal bremsstrahlung emission for 
steep radio spectral indices a > 1. 

It is also possible to explain the hard X-ray and ra- 
dio spectrum by a model where accelerated particles are ex- 
ternally injected in numerous localised regions distributed 
throughout the cluster. In this case, the contribution from 
IC scattering dominates that of non-thermal bremsstrahlung 
in the hard X-ray energy range for steep radio spectral in- 
dices. 

For both the above models a strong assumption was 
made, that the non-thermal X-ray and radio emission from 
Coma is generated from a single electron acceleration mech- 
anism. An advantage of this assumption is that it limits 
significantly the number of free parameters in the model. 
Of course it is possible to have both in situ acceleration 
and external injection of accelerated particles (or global and 
local acceleration) at work in the Coma cluster. We have, 
however, examined the two extremes rather than such a hy- 
brid model. This study has shown that if it is purely in situ 
acceleration that produces the particles necessary for the 
non-thermal emission in the hard X-ray and radio spectral 
bands, then the only explanation for the origin of the hard X- 
ray is bremsstrahlung of supra-thermal electrons; and conse- 
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Figure 6. The energy distribution of electron density for a pure thermal distribution is given by the solid curve, and the dotted curve 
gives best fit to the X-ray (dominated by IC) and radio data for the 'external' source model. The dot-dashed curve gives the combined 
spectrum including the thermal spectrum. For comparison, the best fit in situ model is shown as well (dashed curve). For the 'external' 
model, it is the electrons at ~ 5 X 10 6 keV that are responsible for the hard X-ray (c.f. the in situ model where the hard X-ray is produced 
by bremsstrahlung of the electrons at ~ 50keV). The cluster is assumed to be isothermal with T x = 8.2 keV and an average electron 
density n e ~ 1.23 X 10~ 4 cm~ 3 within the Beppo-SAX field of view . 



quently, to fit both X-ray and radio emission a cluster-wide 
magnetic field of order fj,G (i.e. close to the equipartition 
value) is implied. On the other hand, if the accelerated par- 
ticles are purely provided by numerous, localised external 
sources of injection, then the only possible origin for the 
hard X-rays is inverse Compton scattering of CMB photons 
by relativistic electrons, which implies a cluster-wide mag- 
netic field of 0.1-0.2 fiG. Any independent means of mea- 
suring the global cluster magnetic fields would distinguish 
between the two models. 

At present it is not easy to choose between the models 
of particle acceleration for Coma. Both models: in situ ac- 
celeration (global) and shock wave acceleration (local) have 
their advantages and problems. In favour of the in situ model 
is the similar morphology between the thermal X-ray and 
the radio emission. The radio data have not shown signifi- 
cant structure in the radio intensity in the Coma halo which 
would indicate regions of electron acceleration. Analysis of 
the thermal X-ray emission from the Coma cluster also did 
not show the strong temperature variations expected from 
shocks (Arnaud et al. 2001). On the other hand, several fea- 
tures indicating recent and ongoing mergers are found in the 
Coma cluster (e.g. Donnelly et al. 1999, Arnaud et al. 2001). 
These mergers often produce multiple large scale shocks and 
turbulence capable of both magnetic field amplification and 
in situ electron acceleration (see a recent paper by Ohno et 



al. 2002 on the subject). Whether or not these regions ex- 
tend far enough to allow electrons to fill the whole volume 
of the halo, remains to be discovered. If the origins of the ra- 
dio and X-ray emitting electrons are different, then a much 
wider range of possible models can be accommodated. For 
example, there has been revived interest in the class of sec- 
ondary electron models where relativistic and thermal pro- 
tons in the ICM collide to produce relativistic (secondary) 
electrons (Dennison 1980; Blasi & Colafrancesco 1999; Dolag 
& Enfilin 2000; and Miniati et al. 2001). 
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Figure 7. The ratio of non-thermal bremsstrahlung to inverse-Compton emission at 70 keV derived from the electron energy distributions 
for in situ acceleration at rate aop~ q or the 'external' source model with injection rate Kp~ r > . The q, rj value corresponding to the radio 
halo spectral index is shown as a dotted line. The horizontal dashed line marks when the contributions from non-thermal bremsstrahlung 
and IC are equal. 
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